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Baryon resonances from a novel fat-link fermion action 
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We present first results for masses of positive and negative parity excited baryons in lattice QCD using an 
0(a 2 ) improved gluon action and a Fat Link Irrelevant Clover (FLIC) fermion action in which only the irrelevant 
operators are constructed with fat links. The results are in agreement with earlier calculations of N* resonances 
using improved actions and exhibit a clear mass splitting between the nucleon and its chiral partner, even for 
the Wilson fermion action. The results also indicate a splitting between the lowest J p = | states for the two 
standard nucleon interpolating fields. 



1. INTRODUCTION 

Understanding the dynamics responsible for 
baryon excitations provides valuable insight into 
the forces which confine quarks inside baryons 
and into the nature of QCD in the nonpertur- 
bative regime. One of the long-standing puzzles 
in spectroscopy has been the low mass of the 
first positive parity excitation of the nucleon (the 
J p = | A^*(1440) Roper resonance) compared 
with the lowest lying odd parity excitation. An- 
other challenge for spectroscopy is presented by 
the A 1 / 2- (1405), whose anomalously small mass 
has been interpreted as indicating strong coupled 
channel effects involving T,tt, KN, ■■■ (!]], and 
a weak overlap with a three valence constituent 
quark state. 

In this paper we present the first results of 
excited octet baryon mass simulations using an 
0{a 2 ) improved gluon action and an improved 
Fat Link Irrelevant Clover (FLIC) [|| quark ac- 
tion in which only the irrelevant operators are 
constructed using fat links. Configurations are 
generated on the new Orion computer cluster 
dedicated to lattice gauge theories at the CSSM 
at Adelaide University. After reviewing in Sec- 
tion 2 the main elements of lattice calculations 
of excited hadron masses, we describe in Section 
3 various features of interpolating fields used in 



this analysis. In Section 4 we present results for 
J p = i nucleons and hyperons. Finally, in Sec- 
tion 5 we make concluding remarks and discuss 
some future extensions of this work. 

2. BARYONS ON THE LATTICE 

The history of excited baryons on the lattice 
is quite brief, although recently there has been 
growing interest in finding new techniques to iso- 
late excited baryons, motivated partly by the ex- 
perimental N* program at Jefferson Lab. Previ- 
ous work on excited baryons on the lattice can be 
found in Refs. §-§. 

Following standard notation, we define a two- 
point correlation function for a baryon B as: 

G B (t,p)^J2 (0\Txb(x) xb(P)\0) . « 

where xb is a (positive parity) baryon interpolat- 
ing field, and we have suppressed Dirac indices. 
The choice of interpolating field xb is discussed 
in Section 3 below. For large Euclidean time, the 
correlation function can be written as a sum of the 
lowest energy positive and negative parity contri- 
butions: 
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where a fixed boundary condition in the time 
direction is used to remove backward propagat- 
ing states, and where the overlap of the field Xb 
with positive or negative parity states \B^) is 
parameterized by a coupling strength Xg± , with 



Eb± — y Mg ± + p the energy. The energies 
of the positive and negative parity states are ob- 
tained by taking the trace of Gb with the opera- 
tor T±, where 
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For p = 0, E B ± = M B ± and the operator r± 
projects out the mass, M B ±, of the baryon B ± . 
In this case, positive parity states propagate in 
the 1, 1 and 2, 2 elements of the Dirac matrix of 
Eq. (^), while negative parity states propagate in 
the 3, 3 and 4, 4 elements. 

3. INTERPOLATING FIELDS 

In this analysis we consider two types of inter- 
polating fields which have been used in the litera- 
ture. The notation adopted follows that of Lein- 
weber et al. ||. For the positive parity proton 
we use as interpolating fields: 

Xi + (x) = tabc (ul(x) C7 5 d b (x)) u c (x) , (4) 

and 

X2 + i x ) = £abc(ul(x) C 4(20)75 u c {x) , (5) 

where the fields u, d are evaluated at Euclidean 
space-time point x, C is the charge conjugation 
matrix, a, b and c are color labels, and the super- 
script T denotes the transpose. As pointed out 
by Leinweber ||, because of the Dirac structure 
of the "diquark" in the parentheses in Eq. (|j), the 
field Xi + involves both products of upper x upper 
x upper and lower x lower x upper components 
of spinors for positive parity baryons, so that in 
the nonrelativistic limit Xi + — Further- 
more, since the "diquark" couples to a total spin 
0, one expects an attractive force between the two 
quarks, and hence a lower energy state than for 
a state in which two quarks do not couple to spin 
0. 

The %2 + interpolating field, on the other hand, 
is known to have little overlap with the ground 



state [|3|]9|]. Inspection of the structure of the 
Dirac matrices in Eq. (||) reveals that it involves 
products of upper x lower x lower components 
only for positive parity baryons, so that xf + = 
0(p 2 /E 2 ) vanishes in the nonrelativistic limit. As 
a result of the mixing, the "diquark" term con- 
tains a factor a ■ p, meaning that the quarks no 
longer couple to spin 0, but are in a relative L = 1 
state. One expects therefore that two-point corre- 
lation functions constructed from the interpolat- 
ing field xV~ are dominated by larger mass states 
than those arising from Xi + - 

Interpolating fields for a negative parity pro- 
ton can be constructed by multiplying the posi- 
tive parity fields by 75, x B ~ = 75 X B+ \ which re- 
verses the role of the terms in Eq. (Q) . While the 
masses of negative parity baryons are obtained 
directly from the (positive parity) interpolating 
fields in (|J) and (^) by using the parity projec- 
tors r±, it is instructive nevertheless to exam- 
ine the general properties of the negative parity 
interpolating fields. In contrast to the positive 
parity case, both the interpolating fields Xi~ an d 
X%~ mix upper and lower components, and con- 
sequently both Xi~ and xf~ are ®(p/E)- Phys- 
ically, two nearby J p = | states are observed 
in the nucleon spectrum. In simple quark mod- 
els, the splitting of these two orthogonal states is 
largely attributed to the extent to which scalar di- 
quark configurations compose the wave function 
H . It is reasonable to expect Xi~ to have better 
overlap with scalar diquark dominated states, and 
thus provide a lower effective mass in the large 
Euclidean time regime explored in lattice simula- 
tions. If the effective mass associated with the x\ 
correlator is larger, then this would be evidence 
of significant overlap of x%~ with the higher lying 
Ni~ states. In this event, further analysis di- 
rected at resolving these two states is warranted. 

4. RESULTS 

In this paper we report results of calculations 
of octet excited baryon masses performed on a 
16 3 x 32 lattice at j3 = 4.60 with a lattice spac- 
ing of a = 0.125(2) fm. The analysis is based on 
a preliminary sample of 50 configurations gener- 
ated on the new Orion computer cluster at the 
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Figure 1. Masses of the nucleon and the lowest 
J p = \ excitation. The FLIC and Wilson re- 
sults are from the present analysis, with the NP 
improved clover and DWF |6) results shown 
for comparison. The empirical N and iV*(1535) 
masses are indicated by the asterisks. 



CSSM, Adelaide. For the gauge fields, a mean- 
field improved plaquette plus rectangle action is 
used, while for the quark fields, the FLIC @ 
action is implemented. In the present analysis 
we have imposed fixed boundary conditions in 
the time direction (U t (x,nt) — V x), and pe- 
riodic boundary conditions in spatial directions. 
Although the simulations were performed with 
both n = 4 and 12 fattening sweeps, the im- 
proved gauge fields were found to be smooth af- 
ter only 4 sweeps. Since the results with n = 4 
sweeps exhibit slightly better scaling than those 
with n = 12 ||, we shall focus on the results with 
4 smearing sweeps. The 12 sweep results lead to 
the same conclusions as presented in the follow- 
ing. Further details of the simulations are given 
in Ref. §. 

In Fig. 1 we show the N and N*(h ) masses 
as a function of the squared pseudoscalar meson 
mass, m^. The results of the new simulations 
are indicated by the filled symbols (filled circles 



are FLIC; filled diamonds are Wilson). For com- 
parison, we also show results from earlier sim- 
ulations with domain wall fermions (DWF) j|] 
(open triangles), and a nonperturbatively (NP) 
improved clover action at j3 = 6.2 j?J. The scat- 
ter of the different NP improved results is due to 
different source smearing and volume effects: the 
open squares are obtained by using fuzzed sources 
and local sinks, the open circles use Jacobi smear- 
ing at both the source and sink, while the open 
diamonds, which extend to smaller quark masses, 
are obtained from a larger lattice (32 3 x 64) us- 
ing Jacobi smearing. The empirical masses of the 
nucleon and the lowest \ excitation are indi- 
cated by the asterisks along the ordinate. There 
is excellent agreement between the different im- 
proved actions for the nucleon mass, in particular 
between the FLIC, NP improved clover JtJ and 
DWF § results. On the other hand, the Wil- 
son results lie systematically low in comparison 
to these due to large 0(a) errors in this action 
A similar pattern is repeated for the N*(j ) 
masses. Namely, the FLIC, NP improved clover 
and DWF masses are in agreement with each 
other, while the Wilson results again lie system- 
atically low. A mass splitting of around 400 MeV 
is clearly visible between the N and N* for all 
actions, including the Wilson action, contrary to 
previous claims Q. 

Figure 2 shows the ratio of the masses of the 
A^*(i ) and the nucleon. Once again, there is 
good agreement between the FLIC and DWF 
actions. However, the results for the Wilson 
action lie above the others, as do those for 
the anisotropic D234 action ||. The D234 ac- 
tion has been mean-field improved, and uses an 
anisotropic lattice which is relatively coarse in the 
spatial direction (a w 0.24 fm). This is an indi- 
cation of the need for nonperturbative or fat link 
improvement. 

The mass splitting between the two lightest 
nearby iV*(±~) states (iV*(1535) and AT(1650)) 
can be studied by considering the xi an< i Xi i n_ 
terpolating fields in Eqs.(^) and (0). Recall that 
the "diquarks" in xi and X2 couple differently to 
spin, so that even though the correlation func- 
tions built up from the xi an d Xi fields will be 
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Figure 2. Ratio of the A*(i ) and N masses. 
The FLIC and Wilson results are from the present 
analysis, with results from the DWF 0] and D234 
|Bj actions shown for comparison. The empirical 
iV*(1535)/A mass ratio is denoted by the aster- 
isk. 
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Figure 3. Masses of the 2 

the FLIC action. The positive parity (Ai and A2) 
and negative parity (Aj* and A|) states are con- 
structed from the Xi an d X2 interpolating fields, 
respectively. The empirical masses of lowest two 
\ excitations of the nucleon are indicated by the 
asterisks. 



made up of a mixture of many excited states, they 
will have dominant overlap with different states 
yielding different masses ||. The results, shown 
in Fig. 3 for the FLIC action, indicate that in- 
deed the A* (5 ) corresponding primarily to the 
X2 field (labeled "A|") lies systematically above 
the iV*(i ) associated primarily with the xi ne ld 
("TV*"). 

As has long been known, the positive parity X2 
interpolating field ("A2", which is also sometimes 
denoted by "A'(| )") does not have good over- 
lap with the nucleon ground state ||, and cor- 
responds to a state which lies around 400 MeV 
above the negative parity excitation of xi- There 
is little evidence that this state is the AT* (1440) 
Roper resonance (first i + excitation of the nu- 
cleon). While it is possible that the Roper res- 
onance may have a strong nonlinear dependence 
on the quark mass at < 0.2 GeV 2 , arising 
from pion loop corrections, it is unlikely that this 
behavior would be so dramatically different from 



that of the A* (1535) so as to reverse the level 
ordering obtained from the lattice. A more likely 
explanation is that the X2 interpolating field does 
not have good overlap with either the nucleon or 
the A*(1440), but rather a (combination of) ex- 
cited i state(s). 

Recall that in a constituent quark model in a 
harmonic oscillator basis the mass of the Roper 
is higher than the mass of the lowest P-wave ex- 
citation. The lattice data thus appear to be con- 
sistent with the naive quark model expectation 
at large values of m q . Better overlap with the 
Roper resonance is likely to require more exotic 
interpolating fields. 

5. CONCLUSION 

We have presented the first results for the ex- 
cited baryon spectrum from lattice QCD using 
an 0(a 2 ) improved gauge action and an improved 
Fat Link Irrelevant Clover (FLIC) quark action in 
which only the links of the irrelevant dimension 
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five operators are smeared. The simulations have 
been performed on a 16 3 x 32 lattice at /3 = 4.60, 
providing a lattice spacing of a — 0.125(2) fm. 
The analysis is based on a set of 50 configurations 
generated on the new Orion computer cluster at 
the CSSM, Adelaide. 

Good agreement is obtained between the FLIC 
and other improved actions, such as the nonper- 
turbatively improved clover JtJ and domain wall 
fermion Q actions, for the nucleon and its chiral 
partner, with a mass splitting of ~ 400 MeV. Our 
results for the A*(±~) improve on those using the 
D234 |§ and Wilson actions. Despite strong chiral 
symmetry breaking, the results with the Wilson 
action are still quite reasonable rendering earlier 
conjectures invalid. Using the two standard nu- 
cleon interpolating fields, we also confirm earlier 
observations Q of a mass splitting between the 

two nearby | states. We find no evidence of 
1 4" 

overlap with the ^ Roper resonance. 

We have not attempted to extrapolate the lat- 
tice results to the physical region of light quarks, 
since the nonanalytic behavior of 7V*'s near the 
chiral limit is not as well understood yet as that 
of the nucleon |jl^,0. It is vital that future lat- 
tice N* simulations push closer towards the chi- 
ral limit. On a promising note, our simulations 
with the 4 sweep FLIC action are able to reach 
relatively low quark masses (m q ~ 60-70 MeV) 
already. We have also not addressed the ques- 
tion of to what extent quenching may affect any 
of our results. We naturally expect the effects of 
quark loops to be relatively unimportant at the 
currently large quark masses, although quenching 
may well produce some artifacts as one nears the 
chiral limit. 

For future work, we intend to use variational 
techniques to better resolve individual excited 
states, for instance, those corresponding to the 
Ni and N% fields (using a 2 x 2 correlator ma- 
trix). In order to further explore the origin of 
the Roper resonances, more exotic interpolating 
fields involving higher Fock states, or other non- 
local operators should be investigated. Finally, 
the present N* mass analysis will be extended in 
future to include N — > N* transition form factors 
through the calculation of three-point correlation 



functions. 
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